Hidden quasi one-dimensional superconductivity in Sr2Ru04 
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Using an asymptotically exact weak coupling analysis of a multi-orbital Hubbard model of the 
electronic structure of Sr2Ru04, we show that the interplay between spin and charge fluctuations 
leads unequivocally to triplet pairing which originates in the quasi-one dimensional bands. The 
resulting superconducting state spontaneously breaks time-reversal symmetry and is of the form A ~ 
{Px + ipy) z with sharp gap minima and a d- vector that is only weakly pinned. The superconductor is 
topologically trivial and hence lacks robust chiral Majorana fermion modes along the boundary. The 
absence of topologically protected edge modes could explain the surprising absence of experimentally 
detectable edge currents in this system. 

PACS numbers: 71.10.Hf, 71.10.Fd, 71.27.-|-a, 74.20.Rp, 74.70.Pq 



Introduction - Sr2Ru04 is a layered perovskite mate- 
rial, isostructural to the hole-doped 214 family of cuprate 
superconductors. Below T bOK , it exhibits Fermi liq- 
uid behavior and undergoes a superconducting transition 
at Tc = 1.5K. There is compelling experimental evidence 
which suggests that this superconducting state has odd 
arity[ll III and spontaneously breaks time-reversal [2- 
symmetry. One of the simplest superconducting gap 
functions which meets both of these requirements is the 
chiral p-wave state, A(p) oc {px + ipy) z, a quasi- two- 
dimensional version of superfluid '^He- A [3, Q . 

In its simplest form, this chiral pairing gives rise to a 
topological superconductor: all Bogoliubov quasiparticle 
excitations are gapped in the bulk whereas topologically 
protected chiral Majorana fermion modes exist at the 
edge of the system and in vortex cores . These modes 
are robust against all perturbations, including disorder, 
so long as the BCS pairing gap in the bulk remains finite. 
In addition, spontaneous supercurrents are expected at 
sample edges and domain walls Q . 

However, scanning squid and Hall bar imaging studies 
flol | have revealed that edge currents of the expected mag- 
nitude are not found in Sr2Ru04. Moreover, low temper- 
ature power laws in the electronic specific heatfTlj and 
the nuclear spin relaxation 1/Ti[l3] suggest that this ma- 
terial is not a simple chiral superconductor, which would 
exhibit exponentially activated behavior in both of these 
quantities. The form of the superconducting order pa- 
rameter which accounts for all of the observed phenom- 
ena remains unknown. Resolution of this puzzle could 
come from a careful consideration of the normal state 
properties, which are known with unprecedented detail 

El [3. 

The Fermi surface of Sr2Ru04Consists of 3 sheets, de- 
noted a, P, ^[TT 
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The a and /? sheets are hole and 
electron pockets respectively; they are comprised primar- 
ily of the Ru dxzjdyz orbitals which form quasi one di- 
mensional bands. The 7 sheet is composed mainly of 
the Ru dxy orbital, which forms a quasi two dimensional 
band. A variety of experiments have shown that the sys- 
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FIG. 1: (Color online) Tight-binding Fermi surface for the 
non-interacting Hamiltonian _ffo. Ifybridizations among dif- 
ferent orbitals are neglected in the sold curve and are included 
in the dashed curve. 



tem behaves as a quasi two dimensional Fermi liquid with 
considerable effective mass enhancements [l^ . Therefore, 
it is likely that electron correlations play a significant role 
in influencing the pairing mechanism of this system. 

In this paper, we present a microscopic theory of su- 
perconductivity in Sr2Ru04. Using a simple extension of 
a recently developed asymptotically exact weak-coupling 
analysis of the Hubbard model[l5|, we show that the 
dominant superconducting instability is in the triplet 
channel and occurs on the quasi- ID Fermi surfaces of 
Sr2Ru04. The resulting superconducting state sponta- 
neously breaks time-reversal symmetry. It exhibits node- 
like behavior since it possesses points on the Fermi sur- 
face where the gap is parametrically small. It supports 
Andreev bound states at domain walls and at the edges 
of the system. However, it is a topologically trivial su- 
perconductor without chiral Majorana fermion modes, or 
detectable spontaneous supercurrents at the edges. 

Microscopic model - We consider a simple Hamiltonian 
with three bands derived from the Ru t2g orbitals 



SH (1) 



Here, we introduce vector indices such that a = x, y, and 
z, refer, respectively, to the Ru d^z, dyz and dxy orbitals, 
riiaa is the density of electrons having spin a at position 
i in orbital a and riia = X^o- "-iaa-- The strength of the 
repulsive interaction between two electrons on like (dis- 
tinct) orbitals at the same lattice site is given by U(V). 

Hf) = y^„ y^,r { ~ ~ a] ct cr is the dominant intra- 

orbital kinetic energy and gives rise to three decoupled 
energy bands at the Fermi level as shown in Fig. [1] Here, 
we make use of the following tight-binding parametriza- 
tion of these energies: 



'^l{y)i^) = -2t COS k^f^y 



2t^ cos fct,(a;) 

\{k) = —2t' (cos kx + cos ky) — At" cos k^ cos ky (2) 

where we takc{t,t^ ,t' ,t'\ fi) ^ (1.0, 0.1, 0.8, 0.3, 1.0)[ii 
17|. The quantity 6H represents smaller terms such 
as longer range hopping and spin orbit coupling (SOC) 
which mix the distinct orbitals. It plays a relatively mi- 
nor role in determining the superconducting transition 
temperature. However, SH plays a crucial role in select- 
ing a superconducting state which breaks time-reversal 
symmetry, as will be discussed below. When SH = 0, 
the non-interacting susceptibilities of the normal state 
are separate functions for each orbital: 



d^k /(e^^g+g-) - /(e^^g) 



(3) 



where /(e) is the Fermi function. Since the quasi-two 
dimensional band is almost circular with a radius k'j'^, 
its susceptibility is nearly constant: Xz ~ l/47rt' for q < 
2kJ'^. In contrast, the quasi-lD bands have susceptibili- 
ties that are peaked at = (2fcy'^, tt) and qy = (tt, 2kj'^) 
for the X and y orbitals, respectively. It is the structure 
of Xx and Xy which gives rise to the incommensurate spin 
fluctuations in the material [l^ . 

Since the superconductivity in Sr2Ru04 evolves out of 
a Fermi liquid and Tc ^ i?/, it is reasonable to carry out 
a weak coupling analysis which treats the limit U,V 
where W is the bandwidth. In this limit, superconduc- 
tivity is the only instability of the Fermi liquid, and it 
can be treated in an asymptotically exact manner via a 
two-stage renormalization group analysis [l^. In the first 
stage, high energy modes are perturbatively integrated 
out above an unphysical cutoff, and an effective particle- 
particle interaction in the Cooper channel is derived: 

r,(fc,g,a) = L/ + C/2xa(fc + g)-2l^2^xMfc-g) 
Ttiklc^) = -C/2xa(fc-<7)-2F2^xM^-9) (4) 

where ra(k,q,a) (a = s or t) is the effective interac- 
tion in the singlet (triplet) channel. In the second stage. 
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FIG. 2: (Color online) Pairing eigenvalues as a function of 
V/U for the bandstructure parameters quoted in the text. 
The strongest pairing strengths occur among the quasi- ID 
{a, j3} bands. There is a near degeneracy of the singlet and 
triplet eigenvalues for V = 0, but with V > 0, the quasi-lD 
triplet state is the dominant superconducting configuration. 



the renormalization group flows of these effective interac- 
tions are computed and the superconducting transition 
temperature is related to the energy scale at which an 
effective interaction grows to be of order 1. Following 
this prescription, one obtains 



(5) 



where Aq"'"'' is the most negative eigenvalue of the matrix 



(a, a) 



Vf(k) 



ra(fc,g,a). 



(6) 



with k and q constrained to lie on the Fermi surface of 
band a. The pair wavefunction in the superconducting 
state is proportional to the associated eigenfunction jl5 1 . 



The values of A'"'"-*, obtained by numerical diagonal- 
ization, are presented in Fig. [21 When V = 0, the two 
dimensional z band has its dominant pairing instability 
in the singlet dx2_y2 channel and a substantially lower 
pairing strength in the triplet p-wave channel. By con- 
trast, the pairing tendencies of the x and y bands are 
stronger, and exhibit a close competition between sin- 
glet and triplet pairing [231 . When V > 0, only triplet 
pairing in the quasi-one dimensional bands is enhanced. 
Since the solutions with weaker pairing strengths have ex- 
ponentially smaller transition temperatures in the weak- 
coupling limit, the effect of subdominant orders is negligi- 
ble. Thus, in the asymptotically weak coupling limit, the 
dominant superconducting instability occurs in the quasi- 
one dimensional bands in the spin-triplet channel. For 
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FIG. 3: (Color online) (a) In the absence of any band mixing, 
the triplet state within the xz orbital (red) has kx-w&ve sym- 
metry and has line nodes near ky = 7r/2 (dashed line). The 
condensate on the yz orbital (grey) is related to the one shown 
here by a 90-degree rotation, (b) The chiral state which re- 
sults when small band mixing perturbations are taken into 
account. The relative phase factors on different portions of 
the Fermi surface are shown. 



> 0, triplet pairing in the x and y bands is enhanced 
by virtual charge fluctuations occurring in the z band. 
This conclusion is robust against a rather large range of 
quantitative changes of the band parameters [l9j. 

The triplet pair vifave-function in each band 



X, y 



(7) 



is specified by the complex vector da{k) in spin-space. 
In general, the real and imaginary parts of da are inde- 
pendent real vectors, and the net spin magnetization of 
band a is Ma oc d* x da- However, in weak coupling, 
only "unitary states," i.e. states with Ma = 0, need be 
considered, so da can be expressed, up to a phase, as a 
real "gap function" times a real unit vector 0^: 



da{k) = Aa{k)e"'"n 



(8) 



Figure [3^ shows the sign of Ax{k) on the xz Fermi sur- 
face. In addition to having odd parity, the wave function 
has two point nodes on the Fermi surface near ky = ±7r/2 
and is well approximated by Ax(k) w Ao sinfcj; cosfcy. 
By symmetry, the gap function Ay{k) w Aq sin ky cos k^ 
on the yz Fermi surface. For SH=0, Cl^, (ly, Ox and 9y 
are undetermined. 

Effect of small interactions - We now consider the ef- 
fect of SOC and the hybridization among the different 
orbitals: 

ka 

<^,P k 

where g(k) = —2t" sin k^ sin ky, the second quantity 
above is the SOC, and the angular momentum opera- 
tors are expressed in terms of the totally anti-symmetric 



tensor as i^^ = ieaais- Recent electronic structure cal- 
culations have produced the estimates t" w O.lt and 
77 w 0.1t[il[2^. 

There are several important qualitative effects of in- 
cluding these additional terms in the Hamiltonian: 1) 
Non-zero t" or 77 pins the relative phase, 9^ — 9y, and 
relative orientation, • of the d vectors. 2) Non-zero 
r] defines a preferred ordering direction for Via. 3) The 
nodes on the xz and yz Fermi surfaces are gapped, al- 
though the gap is parametrically small for small SH. To 
understand the role of SH in selecting among the large 
number of possible ordered phases, it is simplest to con- 
sider the Landau free energy J- to low order in powers of 
the order parameter which is a valid approximation near 
Tc. Since any order induced on the z band is slaved to the 
primary order on the x and y bands, we keep explicitly 
only a = X and y, in which case 



J2 [ r\do.\^ + u\da\^ + 7|4 X da\' 



+a^[\dl\' + \dl\']+a,[\d%\' + \dyy\'] 



(10) 



-vi\dx\ \dy\ + V2\dx ■ dy\ 
^ ^ 2 



Ji 



d% ■ dy + dy ■ da 



J2 



d* X dy + dy X d^ 



where the terms in the first line survive the SH — > 
limit, Qj '~ 0{rf' /t^), and Vj and Jj have contributions of 
order (?//t)^ and (t"/t)^. To quadratic order in the order 
parameters and to order (jy/t)^, (i"/t)^, this expression 
is the most general one consistent with symmetry, but for 
the quartic terms, in the interest of simplicity, we have 
assumed that the SOC is weaker than the band-mixing, 
and so have enforced spin rotational symmetry. Since 
\vj \ <^ u, the terms on the third line are not qualitatively 
important. When the action is derived from any form of 
BCS theory, it is possible to show that 7 and Jj > 0. 
Thus, there are two possible phases depending on the 
sign and magnitude of ai: (a) For ai > Min[0,a2], there 
are time-reversal symmetry preserving "B" phase states 
(analogous to the B-phase in '^He) in which fl^ ■ = 



and 6x 



Depending on the sign 



of 02, either fl^. = x or fix ~ y. (b) For ai > Min[0, 02], 
there is a "chiral p+ip" state (analogous to the A phase 
in He3) with fix = ^y = z and 0x — Oy zl 7r/2. All 
other configurations have a higher Free energy and can 
be neglected. 

The parameters aj can be related to differences of sus- 
ceptibilities of the non-interacting system, and so can 
be computed directly from the assumed band-structure. 
For the stated parameters, we find that oi < Min[0,a2], 
so the chiral state is preferred. However, the balance is 
delicate, and this conclusion is not robust against small 
changes to the model. 

Next, we address the fate of the gap nodes when 
SH ^ 0. We have studied the BdG Hamiltonian for 
the chiral state using the gap functions derived in the 
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previous section. Generically, the resulting state is 
nodelcss. However, although the nodes are not topo- 
logically stable, they are parametrically small: where 
a gap node occurs for SH = 0, the induced gap is 
~ Ao [0{t"^/f) + Oljf/t^)] . The energy scale of these 
gap minima is therefore two orders of magnitude smaller 
than the transition temperature. 

Topological •properties and edge currents - Next, we 
consider the topological properties of the system assum- 
ing that ai < Min[0, 02], so the chiral state is preferred. 
The BdG Hamiltonian for the quasiparticle excitations 
in the superconducting state can then be expressed in 
terms of Anderson pseudo-spins as 



* r- 



(11) 



where ^^^j are Nambu spinors, = a, /3 runs over the 
two quasi- ID bands, r are the Pauli matrices, and the 
pseudo-Zceman field is 



S,{k) = Rc[A,(fc)],Im[A,(fc)],e,(fc) - A* 



(12) 



For the chiral p-wave state, the pseudospin has the form 
of a skyrmion in momentum space: it points along the 
— z(-|-z)-direction inside(outsidc) the Fermi surface, and 
on the Fermi surface, it lies in plane, winding by 27r 
around the Fermi surface. The topological properties of 
the chiral state come from the integer skyrmion number 



AC 



1 



47r|<5^(fc)|3 



X dyS„ 



(13) 



where \S^ik)\ = [e^{k) - + \A^{k)\^ . The net 
number of chiral quasiparticle modes at the edge of the 
superconductor is given by the skyrmion number, and so 
long as A/l/ 0, these modes are topologically protected, 
and cannot be localized by backscattering. 

As Afi, is an integer, small changes in the spectrum 
do not affect it. However, it is odd under A^, — > A* 
(i.e. upon transforming px + ipy Px — iPy) or under 
a particle- hole transformation, ~ pL ^ pL ~ e^. As can 
be seen in Fig. 1, hybridization between the two quasi- 
ID bands results in the closed a and j3 Fermi surfaces, 
the former electron-like and the latter hole-like. Con- 
sequently, in a chiral px + ipy state, TVq = —Mp, or in 
otherwords the net skyrmion number is zero so it is not a 
topological superconductor! The chiral edge modes along 
the boundary of the superconductor and along domain 
walls are not protected: in the presence of disorder or 
interactions that scatter a pair from one Fermi surface to 
the other, the counter-propagating edge modes from the 
two bands are localized |19j . 

The existence of edge currents in a chiral p-wave su- 
perconductor follows from symmetry, but whether or not 
it is large enough to detect is still an issue. There is a 



bulk contribution to the edge current which originates 
from the multi-component nature of the order parame- 
ter, and one from the chiral quasiparticle modes near the 
edge of the system. We will address the question of how 
the highly non-circular band structure of the quasi ID 
bands affects the bulk contribution in a future publica- 
tion [lit ■ Since the chiral quasiparticle modes travel with 
velocity v = Vf{A/Ef), they would, by themselves, make 
a readily detectable contribution 0{A/Ef) to the edge 
current, were they not localized. 

Discussion - It follows from general arguments plj 
that near Tc, superconductivity can arise either in the 
{a, 13} pockets or in the 7 pocket; below Tc, supercon- 
ductivity is induced in the subdominant Fermi surfaces 
via a proximity effect 



(14) 



with J' <C J7(2lJ. Due to the weakness of this proximity 
effect, it can be expected that for a range of tempera- 
tures J'{A/Ef) <C r < Tc, superconductivity is present 
essentially only on the dominant portions of the Fermi 
surfaces. Nonetheless, in the present case, this coupling 
produces a small gap on the xy Fermi surface, which then 
adds to the skyrmion number; while the energy scales in- 
volved are likely too small to affect the results of any 
practical experiment, ultimately the proximity effect re- 
stores the topological character of the chiral supercon- 
ducting state. 

As we have shown, it is unavoidable, given the band- 
structure of Sr2Ru04, and assuming the interactions are 
weak, that the superconductivity arises primarily on the 
quasi-lD bands. Within this framework, there is a natu- 
ral explanation for the surprising absence of superconduc- 
tivity in the closely related bilayer compound Sr3Ru207; 
the bilayer splitting in Sr3Ru207 primarily affects the 
quasi-lD bands, leaving the 2D dxy band essentially un- 
changed. If pairing occurred primarily on the quasi-2d 
band, both compounds ought to be similarly supercon- 
ducting. (An analysis of superconductivity in the bilayer 
system will be presented in a separate publication figj ) . 

Admittedly, we cannot rule out the possibility that 
strong coupling effects will change this conclusion. How- 
ever, the absence of experimentally detectable edge cur- 
rents is difficult to reconcile with a state which primarily 
involves the 2D bands. Further analysis of the experi- 
mental evidence concerning this issue will be postponed 
to a future paper [l^. We will also discuss the nature 
of the collective modes which refiect the near indepen- 
dence of the superconducing state on each of the quasi-ld 
bands, especially near the transition temperature. 
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